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Low-temperature scanning electron microscopy (LTSEM) has been used to image the supercurrent distribu- 
tion in ramp-type Josephson junctions between Nb and either the electron-doped cuprate Nd2-,vCe.vCu04_y or 
the hole-doped cuprate YBa2Cu307. For zigzag-shaped devices in the short junction limit the critical current 
is strongly suppressed at zero applied magnetic field. The LTSEM images show, that this is due to the Joseph- 
son current counterflow in neighboring and it facets, which is induced by the d^2_^.2 order parameter in the 
cuprates. Thus, LTSEM provides imaging of the sign change of the superconducting order parameter, which 
can also be applied to other types of Josephson junctions. 
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One of the most controversial topics on high-T^ cuprate su- 
perconductors has been the determination of their order pa- 
rameter symmetry (OPS). A myriad of experiments have been 
performed, indicating a predominant d^i_^i OPS, which im- 
plies important consequences for the microscopic mechanism 
of Cooper pairing in these materials. Obviously, it was quite 
difficult to identify an unambiguous experiment for the de- 
termination of the cuprate OPS. Among the most convinc- 
ing experiments is the observation of half-integer magnetic 
flux quanta in tricrystal grain boundary Josephson junctions 
(JJs) by scanning SQUID microscopy yj. These experi- 
ments, and related integral measurements of critical current 
Ic vs applied magnetic field B, rely on the difference 7i of the 
phase of the order parameter between orthogonal directions 
in (^jc,A;y)-space, which can be detected by interferometer- 
type configurations, such as corner junctions 0, tricrystal 
rings and long JJs (3] H S H ID, and dc 71 SQUIDs 
in [To] [m [121, or by the angular dependence of Ic in biepi- 
taxial JJs fTT\. High-quality hybrid ramp-type JJs, combin- 
ing an s-wave superconductor (Nb) with either the hole-doped 
cuprate YBa2Cu307_s (YBCO) (II [B] [161 or the electron- 
doped cuprate Nd2-.vCevCu04_y (NCCO) |T7l have also been 
realized. Arranging such JJs in a zigzag geometry with the 
facets oriented along the a- and fe-axis of the cuprate, one ob- 
tains alternating facets of and n JJs lITSl [171 . 7t JJs [18 | 
have negative Ic, i. e., js = — y'c sin(|) — 7cSin((|) + 7t), instead 
of is — 7'cSin(|), where js is the supercurrent density; jc > 
is the maximum supercurrent density, and (|) is the Joseph- 
son phase. Realizations include JJs with magnetic baiTiers 
|[T9l|20l|2TJ|22l|23l, geometric constrictions in d-wave super- 
conductors 1241 , Nb JJs with a mesoscopic Au control chan- 
nel li25l , Al JJs with a controllable quantum dot in a In As 
nanowire li26l . and gate-controlled carbon nanotube JJs li27l . 



JJs containing both, 0- and 7t-parts have also been realized 
using ferromagnetic barriers ESl l29l [30l or current injectors 

EH- 

A striking property of s-d-wave zigzag JJs in the /ong JJ 
limit (facet length a~4A,j) is the spontaneous generation of 
magnetic flux ±<I>o/2, i.e. a semifluxon at each corner of 
the zigzag (<I>o = /i/2e is the magnetic flux quantum and 

— 1/2 

'^j°^Jc the Josephson penetration depth). The presence of 
semifluxons in such devices was demonstrated |16| by scan- 
ning SQUID microscopy. In the short JJ limit (neglecting self- 
field effects), for a JJ with facets, the supercurrent density 
in the n''' facet can be described as 1 15| 

;,(i)-(-l)";,(i)sin{(^o + (27t<I>//<I>oA^fl)-i} . (1) 

Here, x is the coordinate along the zigzag (with i = at the JJ 
edge), and <!>/ is the magnetic flux per facet. As the prefactor 
(— 1)" changes sign at every corner of the zigzag, as a direct 
consequence of the c/-wave OPS, Ic{B) is not Fraunhofer-like; 
instead, it has main maxima (/"'"^ — {yiijNjcha for jc{x) = 
const.) at finite field, corresponding to <Pf = ±<I>o/2 for even 
A^, with junction area h ■ a per facet. According to Eq. ([T]l, at 
such <I>^, js{x) — jc\ sm'^^al in each facet. Ic{B) atB — Q has a 
minimum (for even A^) or a small local maximum (for odd A^). 
In the case of homogeneous jc{x), absence of self-field effects 
and even A^ one expects /c(0) — 0, due to a current distribution 
js ^ { — l)"jc, and current reversal at each corner of the zigzag 
results in a quite unusual Ic{B) dependence 031 [TTl , which 
provides strong (indirect) evidence of the Josephson current 
counterflow as a direct consequence of the sign change in the 
d-wa\e order parameter. 

In this Letter we show that low-temperature scanning elec- 
tron microscopy (LTSEM) allows imaging of the supercur- 
rent distribution in YBCO-Nb and NCCO-Nb JJs, and we 
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demonstrate Josephson current counterflow in 0- and 7t-facets 
in zigzag shaped cuprate/Nb JJs at B = 0. 

We investigated hybrid ramp-type JJs with 150 nm thick 
[001] YBCO or optimally doped (;c=0.15) NCCO bottom 
electrodes, grown epitaxially on [001] SrTiOa (STO) single- 
crystal substrates and covered by an STO film with thick- 
ness 100 nm and 35 nm, respectively. After milling a shallow 
ramp (15°-20°) into the bilayers, an epitaxial YBCO (6nm) 
or NCCO (12 nm) interlayer was grown, followed by in-situ 
deposition of a Au barrier layer of thickness Jau, and a Nb 
layer (140-160 nm) as a counter electrode 1 14l [T5l[T7l . In to- 
tal, we investigated four chips with identical layout. Three 
chips contained YBCO-Nb JJs with t/Au=14 nm (chip Yl) and 
12 nm (Y2 and Y3) in order to investigate samples with differ- 
ent jc, i.e. different Xj. The chip N with the NCCO-Nb JJs had 
t/Au=12nm. Below, we show data from chips N, Yl and Y2 
for zigzag JJs with N — % and fl=25 /jm (chip Yl) or A^=10 and 
fl=40 /jm (chip N) and for reference single facet JJs (fl=50 /jm), 
oriented along the a,b axis of the cuprate film (chips N and 
Y2). The conversion from B (normal to the substrate plane) 
to magnetic flux <I> in the JJ was done by comparing the mea- 
sured Ic{B) with /,.(4>) calculated from Eq. (3) in Ref. |15I . 
Considering the idle region (overlap of the Nb electrode on top 
of the cuprate), we can only give a rough estimate on an upper 
limit for the normalized JJ length Na /Xj~2, i.e. all devices are 
expected to be in the short JJ limit. Regarding further electric 
transport properties of our samples, see Refs. lfT4l[T5l[T7l[32l . 

For imaging by LTSEM, the sample was mounted on a He 
cryostage and operated at a temperature T «5-6K. The lo- 
cal perturbation by the focused electron beam (e-beam) cen- 
tered at the position (xo,yo) on the sample surface in the {x,y) 
plane induces an increase in temperature 5T {x — xo,y — yo) on 
a lateral length scale of «l-3/jm, which determines the spa- 
tial resolution of this imaging technique. The maximum local 
increase in temperature Ar is typically < 1 K, and can be ad- 
justed by the e-beam voltage Vt, and beam current 4 ll33l[34l . 
For the LTSEM images shown below Vi,=10 kV and 4=50 pA- 
1 nA. dT results in a local reduction of jc{T) and a concomi- 
tant change of the overall 4 of the JJ. It has been shown the- 
oretically 135] |36| and experimentally [|32l |38| that this ef- 
fect can be used to image the spatial distribution of the super- 
current density js{x) (at / = 4, convoluted with the 5T pro- 
file) along a short JJ by recording the beam-induced change 
blc{x) o= js{x) of the overall critical current as a function of 
the beam coordinate i, during scanning along the JJ. For sim- 
plicity, rather than detecting 54 , we current bias the JJ slightly 
above 4 (typically at a voltage V of a few /jV) and detect the 
beam-induced voltage change bV [34]. Assuming a constant 
differential resistance Rii yields bV{x) = — 7^^/54 °= js{x)- 
To improve the signal-to-noise ratio, we modulate the e-beam 
at 5 kHz (6.6 kHz) and lock-in detect the voltage response 
from the YBCO(NCCO)-Nb JJs. 

In order to characterize the quality of our devices and to 
demonstrate imaging of the current distribution by LTSEM, 
we first present results from the YBCO-Nb and NCCO-Nb 
single facet (A^=l, 0=50 /jm) reference JJs. The inset in 



Fig. JTJa) shows an SEM image of the NCCO-Nb JJ. Fig- 
ure [TFa) shows normalized critical current 4/^™"" vs applied 
magnetic flux <I> = A^<t>/. Fraunhofer-like 4 oscillations are 
clearly visible, although deviations from the ideal characteris- 
tic (dashed line) are obviously present. Those deviations are 
probably mainly due to the finite voltage criterion for the de- 
tection of 4, however also indicate inhomogeneities in jc{x). 

Figures [l|b){T];g) show LTSEM images 5y(xo,yo) for both 
reference JJs (left: YBCO-Nb; right: NCCO-Nb) taken at dif- 
ferent values of <I> as indicated in Fig. [Tf a). At <I> = [main 
4 maximum; graphs (b) and (c)] the voltage signals at yo = 
(i-axis) are positive along the entire length of both JJs. At 
<I> = 0, for N = 1 one finds from Eq. ([T]) that the supercur- 
rent density at 4 is js{x) = jc{x), and hence 5y(i) °^ jc{x), 
i. e., the variation in SV{x) along the JJ directly yields the 
variation of jc{x). The observed 5y(i) clearly indicates 4 
inhomogeneities along the JJs, which are most likely due to 
variations in the quality of the interface and in the thickness 
of the Au barrier layer. For the YBCO-Nb JJ, we find a max- 
imum variation in jc{x) of ±15%. For the NCCO-Nb JJ we 
observe a steplike decrease of jc{x) at i « 35^m by ^30 %. 

The second row of LTSEM images [graphs (d) and (e)] are 
taken at the first side maximum in 4('f'), i- e. at <I> = |<I>o for 
which one expects a sinusoidal variation of the supercurrent 
density js{x) = 4(-^) sin (37ti/fl) with 3/2 wavelengths. This 
behavior is well confirmed by the LTSEM images. The low- 
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FIG. 1: (color online). Single facet NCCO- and YBCO-Nb JJs: (a) 
Normalized critical current Ic/I'"'" vs magnetic flux <I'/'I>o. Labels 
(b)-(g) indicate working points for LTSEM images below. Inset: 
image of NCCO-Nb JJ; dashed frame indicates size and position of 
LTSEM images (b) - (g). Numbers indicate full range |5Vmflx| (in 
^V) of the scale bar (symmetric about 8V = 0). (h),(i): linescans 
5V(x) at <5/<I>o = 5/2 along the JJs, respectively, from images (f) 
and (g), and calculated current density js{x) / jc (solid black lines). 
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est row of LTSEM images [graphs (f) and (g)] for <I> = §<t>o, 
i. e. taken at the second side maximum in again clearly 
shows the expected oscillation with 5/2 wavelengths. The 
graphs (h) and (i) in Fig. [T] show linescans taken from the 
corresponding LTSEM images (f) and (g), together with the 
calculated normalized current density distribution js{x)/jc, 
which was convoluted with a Gaussian beam-induced tem- 
perature profile with a ~ 2.5/jm. The excellent 
agreement between the measured voltage signals and calcu- 
lated current distribution clearly demonstrates that we indeed 
image the supercurrent density distribution along the JJs. 

In the following, we present results on the zigzag JJs [c.f. a 
schematic view in the inset of Fig.j2|a)], starting with the 
YBCO-Nb JJ (N=8, a=25 /jm); Fig. |2[b) shows an SEM im- 
age of this device. Figure |2| a) shows Ic{B) measured on the 
LTSEM cryostage at rw6K (dots) and in a liquid He cryo- 
stat at r=4.2K (solid line). As expected for an array of 0-Jt 
facets, Ic{B) shows main maxima at finite field (B,„ax=i-i jwT) 
and only a small central maximum at B=Q. Because of the 
higher temperature of the LTSEM cryostage, the values are 
reduced, as compared to the 4.2 K data and the oscillations 
are washed out. Nevertheless, almost all maxima and minima 
in Ic{B) still show up at T w 6K. 

For each point of the Ic{B) dependence at 6 K in Fig. 2ta) 
LTSEM images were recorded. Figures |2|c), |2|e) and 2|g) 
[left row] show images taken at three values of B [as labeled 
in graph (a)], namely at the small maximum in Ic{B) at B=Q 
(c), at the main maximum in Ic {B) (e), and at the next side 
maximum in Ic{B) (g). To the right of each LTSEM image, 
we show the corresponding image js{xo,yo) of the supercur- 



rent density distribution (normalized to a spatially homoge- 
neous jc) which was calculated as follows: The ID distri- 
bution js{x) along a zigzag line in the x-y plane was calcu- 
lated numerically from Eq. ([T]i, and all the points {x,y) outside 
the zigzag line were set to y'j = 0. The resulting 2D js{x,y) 
distribution was then convoluted with a Gaussian profile, 
i. e., Mxo,yo) = I^Z ^■^•(^'^) ^M-ryio^jdxdy, 
with = (x — xq)^ + {y — yo)^ and a = 2.5 /jm, and plotted 
in Figs. |2|d), |2|f), |2|h)- The calculated images are in good 
qualitative agreement with the LTSEM images. As the main 
result, Fig.|2|c) clearly shows the alternating sign of supercur- 
rent flow across neighboring facets at B=Q. Thus, the LTSEM 
image provides a direct proof of the existence of and K facets 
in the zigzag JJ, due to the sign change of the order parame- 
ter in the t/-wave cuprate superconductor YBCO. In contrast. 
Fig. |2je) taken at the main maximum in Ic{B), shows only 
positive voltage signals which are largest inside the facets and 
which tend to zero at the corners. This is in qualitative agree- 
ment with js{x) o= |sin7ti/a| as expected for a homogeneous 
zigzag JJ with jc=const. Quantitative differences as observed 
by the LTSEM voltage signals can most likely be attributed 
to jc inhomogeneities along the zigzag JJ, as such inhomo- 
geneities have already been observed for the YBCO-Nb ref- 
erence JJ [c. f. Fig.[TJb)]. The LTSEM image recorded at the 
next side maximum in Ic{B) [Fig.|2|g)] shows a polarity of the 
voltage signals (positive outside and negative in the center) 
which is reminiscent of the behavior of the reference JJs also 
biased at the first side maximum in IciB) [c. f. Figs.[T]:d) and 
[TJe)]- Again, this is in qualitative agreement with the calcu- 
lated js{x) for the zigzag JJ with homogeneous jc distribution 
[c.f.Fig.|2|h)]. 




FIG. 2: (color online). YBCO-Nb zigzag JJ: (a) Ic{B) patterns; inset: 
sketch of zigzag-shaped ramp JJ. (b) Surface image and (c), (e), (g) 
corresponding LTSEM images (Ij, = 50 pA) taken at different values 
for B as indicated in (a); (d), (f), (h) show corresponding calculated 
images of current distribution along the zigzag. 




FIG. 3: (color online). NCCO-Nb zigzag JJ: (a) lc(B) pattern; (b) 
surface image; (c), (e), (g) LTSEM voltage images (/;, = 1 nA) taken 
at different values for B as indicated in (a), (d), (f), (h): correspond- 
ing calculated images of current distribution along the TA%Td%. 
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Finally, we demonstrate that similar results were obtained 
by imaging the current distribution in the NCCO-Nb zigzag 
JJ (A^=10, a=40fjmy, c. f. the SEM image in Fig. [Sjb). Fig- 
ure [3j a) shows Ic{B) measured on the LTSEM cryostage at 
rw5 K, which was almost identical to Ic{B) measured in liq- 
uid He at 4.2 K. As for the YBCO-Nb zigzag JJ, Ic{B) has a 
small central Ic maximum and main maxima at finite field. 
Figures|3jc),|3je) and|3jg) show LTSEM images taken at three 
values of B [as labeled in graph (a)], namely at the small cen- 
tral maximum in Ic{B) at B = (c), at the "dip" in Ic{B) close 
to B=0 (e) and at the main maximum in Ic{B) (g). As in Fig.|2] 
the corresponding calculated images (d), (f), (h) of js{xo,yo) 
(with o=2.5 fjm) are in qualitative agreement with the LTSEM 
images. Again, Fig.[3]^c) clearly shows the alternating sign of 
supercurrent flow across neighboring facets at B=0. This pat- 
tern remains almost unchanged in a small applied field [bias 
point "e" in (a)] as shown in Fig.|3je). Here the polarity of the 
LTSEM voltage signals for the two facets at the right edge of 
the JJ changed. Probably due to the jc inhomogeneity along 
the entire JJ this state results in an even lower value of as 
compared to the Ic value at B=Q. At the main maximum, the 
LTSEM image in Fig.[3]^g) again shows only positive voltage 
signals, as expected, and as discussed above. 

In conclusion, we have shown that low-temperature scan- 
ning electron microscopy allows imaging of the supercurrent 
distribution in cuprate-Nb hybrid ramp-type Josephson junc- 
tions. LTSEM images recorded at B = show Josephson cur- 
rent counterflow. This gives direct evidence of the presence of 
alternating and n facets in YBCO-Nb and NCCO-Nb zigzag 
junctions, which is due to the sign change of the cZ-wave or- 
der parameter in the cuprate superconductors involved in this 
study. We note, that the same technique can also be applied to 
other systems which produce 0-7t Josephson junctions, e.g. JJs 
with a ferromagnetic barrier. Furthermore, this technique may 
also be applied to investigate the order parameter symmetry 
in less studied superconducting materials, if they can be com- 
bined with an i-wave superconductor to form hybrid Joseph- 
son junctions. 

This work was supported by flie DFG (K1930/1 1), the FOM, 
the NWO and by flie ESF programme NES. 
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